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We investigate properties of exotic pairing states in a three-dimensional Fermi gas with three- 
dimensional spin-orbit coupling and an effective Zeeman field. The interplay of spin-orbit coupling, 
effective Zeeman field and pairing can lead to first-order phase transitions between different phases, 
and to interesting nodal superfluid states with gapless surfaces in the momentum space. We then 
demonstrate that pairing states with zero center-of-mass momentum are unstable against finite 
center-of-mass momentum Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) states, with the center-of-mass 
momentum of the pairs opposite to the direction of the effective Zeeman field. Unlike conventional 
FFLO states, these FFLO states are induced by the coexistence of spin-orbit coupling and Fermi 
surface deformation, and have intriguing features like first-order transitions between different FFLO 
states, nodal FFLO states with gapless surfaces in momentum space, and exotic fully gapped FFLO 
states. With the recent theoretical proposals for realizing three-dimensional spin-orbit coupling in 
ultracold atom gases, our work is helpful for the future experimental studies, and provides valuable 
information for the general understanding of pairing physics in spin-orbit coupled fermionic systems. 
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I. INTRODUCTION 

Ultracold atom gases have been considered as ideal 
platforms for the quantum simulation of interesting mod- 
els in a variety of different physical contexts, ranging 
from condensed matter physics, to nulcear physics yj- 
With powerful tools like Feshbach resonance, optical 
lattice etc., strongly interacting strongly correlated mod- 
els can be studied experimentally using dilute ultracold 
atom gases. Of particular importance is the possibility 
to simulate models that have no counterparts in other 
systems. In the past decade or so, interesting phenom- 
ena like BCS-BEC crossover, superfluid to Mott insulator 
transition in a Bose-Hubbard model etc., which do not 
occur naturally in condensed matter systems, have been 
investigated experimentally using ultracold atom gases 

mm- 

Recently, by engineering the atom-laser coupling, syn- 
thetic spin-orbit coupling (SOC) has been realized exper- 
imentally in ultracold atoms |^]7]. As SOC is considered 
to play a key role in systems with interesting properties 
like quantum spin Hall effect or topological phases [8j, 
the experimental realization of SOC in ultracold atoms 
has greatly extended the possibility of quantum simu- 
lation in these systems. A great amount of theoretical 
and experimental efforts have since been dedicated to 
the characterization of rich physics in spin-orbit coupled 
atom gases [5ti30j. Most of these studies have focused on 
the Rashba or Dresselhaus SOC, which also exist in solid 
state systems; or on the NIST SOC, which can be readily 
implemented experimentally [1J|3T]. More exotic forms of 
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SOC, such as a three-dimensional (3D) SOC, which can- 
not be found in natural condensed matter systems, have 
also attracted some attention recently in bosonic systems 
[32"- 38]. With the recent proposals of implementing 3D 
SOC in ultracold atom gases [39 -42 , it is hopeful that in- 
teresting physics like 3D topological insulators |S], Weyl 
semi-metals [43^, or pairing physics in a Fermi gas etc. 
|30| . could be investigated in these systems. 

Indeed, a particularly interesting problem that has 
been under intensive theoretical study lately is the pair- 
ing superfluidity in an attractively interacting Fermi 
gas with SOC and effective Zeeman fields. In a two- 
dimensional (2D) Fermi gas with Rashba SOC and ax- 
ial Zeeman field, it has been shown that a topologi- 
cal superfiuid state can be stabilized, which supports 
topologically protected edge states and Majorana zero 
modes at the core of vortex excitations |5J [TUJ [ST] Hi] . 
While for a 2D Fermi gas with the experimentally avail- 
able NIST SOC or for a 3D Fermi gas with Rashba 
SOC, nodal superfiuid states with gapless excitations ex- 
ist |ia Hg [ig 1211 1271 HH]. Furthermore, it has been 
suggested that for a 2D Fermi gas under NIST SOC 
and cross Zeeman fields, finite center-of-mass momentum 
Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) states become 
the ground state of the system [28,. These FFLO states 
are different from the conventional FFLO pairing states 
in a polarized Fermi gas in that they are driven by SOC- 
induced spin mixing and transverse-field-induced Fermi 
surface deformation |23| [28H30| . A natural question that 
follows is whether such an exotic pairing state is a unique 
feature of systems with SOC and Fermi surface asymme- 
try. It is therefore desirable to study pairing states in 
related systems. 

In this work, we study the pairing states in a 3D at- 
tractively interacting Fermi gas with 3D SOC and an 
effective axial Zeeman field. As previous studies with 
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3D SOC have mostly focused on bosonic systems, our 
work provides a starting point for the understanding of 
the effects of 3D SOC in a fermionic system. Addition- 
ally, the Fermi surfaces of such a system is asymmetric 
along the direction of the effective Zeeman field. Hence 
it is an ideal system for the study of pairing states under 
SOC and Fermi surface asymmetry. To get a qualitative 
understanding of the pairing physics, we focus on the 
zero-temperature phase diagram on the mean-field level. 

For the zero center-of-mass momentum pairing states, 
we find that exotic nodal superfluid states exist with ei- 
ther two or four closed gapless surfaces in momentum 
space. As the anisotropy of the 3D SOC along the axial 
direction increases, these gapless surfaces shrink gradu- 
ally, and will eventually collapse into gapless points as the 
3D SOC is reduced to the Rashba SOC. Therefore, these 
nodal superfluid states can be seen as the counterparts of 
the nodal superfluid states in a 3D Fermi gas with Rashba 
SOC and an effective Zeeman field [16J. We also find 
that the interplay between SOC, effective Zeeman field 
and pairing in the system can lead to first-order phase 
boundaries between different superfluid states, which im- 
plies the existence of phase separated state in a uniform 
gas and spatial phase separation in a trapped gas. These 
findings provide a valuable starting point for understand- 
ing finite center-of-mass momentum pairing states. 

When the finite center-of-mass momentum pairing is 
taken into account, the phase diagram is qualitatively 
modified. Similar to the effects of a transverse field in the 
NIST scheme, the effective axial Zeeman field in combi- 
nation with the 3D SOC makes the Fermi surface asym- 
metric. As a result, pairing states with zero center-of- 
mass momentum become unstable against FFLO states, 
with the center-of-mass momentum of the pairs deter- 
mined by the Fermi surface asymmetry [5S] . These FFLO 
states are characteristic of pairing states in the presence 
of spin-orbit coupling and Fermi surface deformation, 
which induce a shift of the local minima in the ther- 
modynamic potential landscape onto the plane of finite 
center-of-mass momentum. Consequently, these FFLO 
states retain many features of their zero center-of-mass 
momentum counterparts, e.g. the existence of first-order 
boundaries between different FFLO states, the stabiliza- 
tion of nodal FFLO states with gapless surfaces in mo- 
mentum space, and the presence of a continuous transi- 
tion between the gapless and a fully gapped FFLO state. 
As similar scenarios have also been reported in Fermi 
gases with NIST SOC in different dimensions [Ml l29] . 
we expect that these should be the unique features for 
pairing states under SOC and Fermi surface asymmetry, 
and should be independent of the concrete form of SOC. 

The paper is organized as follows: in Sec. I, we present 
our mean-field formalism; in Sec II, we first discuss the 
phase diagram when only zero center-of-mass momentum 
pairing states are considered; we take finite center-of- 
mass pairing states into account in Sec. Ill, where we also 
discuss the general physical picture of pairing under SOC 
and Fermi surface asymmetry. Finally, we summarize in 
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II. MODEL 

We consider a 3D uniform two component Fermi gas 
with 3D spin-orbit coupling crj^ • k_L -I- ^(Jzkz, where ai 
{i — x,y,z) are the Pauli matrices, and the parameter 
7 depicts the anisotropy of the 3D SOC. The system is 
further subject to an effective Zeeman field along the z 
direction. The Hamiltonian of the system can be written 
as 

Ho = Xl('^k-M)(alc«k + 6k&k)- ^Xl("k"k-&k&k) 

k k 

+ ^ ^ ajj.6^i^6_k'ak' + 7 ^ afc cos 6'k(aj,.ak - &k&k) 

k k 

+ ^afcsin0k(e-**-a],5k + e"^''&],ak), (1) 

k 

where k = fc(cos0k,sin0kCOS(/)k,sin0kSin0k), V is the 
quantization volume in 3D, ek = h^k"^ /{2m), Ok (flk) and 
6k (^k) represent the annihilation (creation) operators for 
atoms of different spin species, /i is the chemical poten- 
tial, h is the effective Zeeman field, a is the spin-orbit 
coupling strength. We have assumed s-wave contact in- 
teraction between atoms of different spin species, where 
the bare interaction rate U should be renormalized fol- 
lowing the standard relation 

l/C/ = l/C/p-l/y^l/2ek. (2) 

k 

The physical interaction rate Up is defined as Up — 
ATrh'^as/m, with Qs the s-wave scattering length between 
the two spin species, which can be tuned via the Fesh- 
bach resonance technique. We note that while the SOC 
becomes isotropic when the anisotropy parameter 7 = 1, 
it reduces to the Rashba SOC for 7 = 0. 

The Hamiltonian ([T|) can be diagonalized under the 
mean-field approximation, for which we define the Fulde- 
Ferrell (FF) pairing order parameter [45 : Aq = 
J7^j^(6_k+Q/2'2k+Q/2)- The resulting effective Hamilto- 
nian can be arranged into a matrix form under the basis 

'0k+Q/2 = (ak+Q/2,&k+Q/2,alk+Q/2'^-k+Q/2)^ 

HcS = ^Xl'^I+Q/2^^^^k+Q/2 +Xl(^k+Q/2 -M) 
k k 

- ^IAqP, (3) 

with 

Q _ f Ho{k + Q/2) zAct^y \ 

" -^Aqa, -i/*(-k + Q/2) j ■ 

Here, -ffo(k) — (ek~M)-^~(t~o;7^2)'''z+Q^'^J-'k_L, where / 
is the identity matrix. Diagonalizing the effective Hamil- 
tonian ([3]), we may then evaluate the zero-temperature 



thermodynamic potential 




- Z^IAqP- (5) 

Here the quasi-particle (hole) dispersion ^ — ±) 

are the eigenvalues of the matrix in Eq. (|3| . With- 
out loss of generality, we assume h > 0, Aq = A7and Aq 
to be real throughout the work. 

Typically, the ground state of the system can be found 
by solving the gap equation di^/dA = and the number 
equation n = ^dfl/dfj, simultaneously, where n is the 
total particle number density. However, the competition 
between pairing and the effective Zeeman field leads to 
double-well structure in the thermodynamic potential for 
both the Q = and the finite Q cases. Therefore, one 
must explicitly take the possibility of phase separation 
into account to get the correct ground state for a uni- 
form gas As we will show later, this is also true 
for calculations in a canonical ensemble, when one fixes 
the total particle number and minimizes the Helmholtz 
free energy. Alternatively, we fix the chemical potential fj, 
and look for the global minimum of the thermodynamic 
potential ([s]) as a function of the pairing order parame- 
ter Aq and center-of-mass momentum Q. The resulting 
state corresponds to the local ground state in a trapped 
gas under the local density approximation (LDA). The 
total particle number in the trap can then be evaluated 
by integrating the local number densities from the trap 
center to its edge, i.e. by integrating the number equa- 
tion with decreasing chemical potentials. One advantage 
of our algorithm is that the resulting phase diagram can 
be directly mapped to the phase structure in slow- varying 
potential traps under LDA. 

III. PAIRING STATES WITH ZERO 
CENTER-OF-MASS MOMENTUM 



3 



k 




-1.5 -1.0 -0.5 0.5 1.0 1.5 ' -1.5 -1.0 -0.5 0.5 1.0 1.5 



^7 




-1.5^ -IS 



Figure 1: Gap less contours in the kx — kz plane for 7 = 0.5 
and (fchOs)"^ = —1, with the parameters: (a) akh/h — 0.85, 
fi/h = 0.8; (b) akh/h = 1, fi/h = 0.8; (c) akh/h = 1.05, 
fi/h = 0.8; (d) akh/h = 1.25, fi/h = 0.8; (e) akh/h = 1.25, 
fi/h = 0.35; (f) akh/h = 1.25, fi/h = 0.2. Note the gapless 
contours have axial symmetry around z-axis. The effective 
Zeeman field h is taken to be the unit of energy, while the 
unit of momentum kh is defined through f?k\/2m = h. 



We start by analyzing the pairing states with zero 
center-of-mass momentum, which will provide us with 
valuable understandings for the FFLO pairing states that 
we will discuss in the following section. For Q = 0, an 
analytical expression for the quasi-particle (hole) disper- 
sions is typically not available, except for those along the 
kz axis, which can be written as 



Et=o..± = -^[2 ± ^(^k - /i + Xo^lKY + A2] . (6) 

Apparently, two of the branches _ can cross zero, 
leading to nodal superfluid states with gapless excita- 
tions. As the dispersion spectra have the symmetry 
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the gapless points on the kz axis must 
be symmetric with respect to the origin. Furthermore, we 
find that the gapless points in momentum space typically 



form closed surfaces that are rotationally symmetric with 
respect to the kz axis. Typical gapless contours in the 
k'x — kz plane are shown in Fig. ([T| , where we demonstrate 
the evolution of the gapless surfaces as the parameters 
are tuned. It is clear that there can be two or four sep- 
arate gapless surfaces in momentum space, for which we 
may define different nodal superfluid states. The num- 
ber of these closed surfaces in momentum space can be 
deduced by counting the number of gapless points on the 
kz axis. From Eq. Q, we know that there are eight 
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Figure 2: Phase diagram obtained in the fj, — a plane with (a) 
(fchOs)"^ — 0, and (b) {khaB)~^ = —1, respectively. The bold 
solid curves represent the first-order boundaries, and the thin 
dashed curves represent continuous phase boundaries. For 
both figures, we have set 7=1. 
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no gapless points otherwise. These correspond to nodal 
superfluid state with four closed gapless surfaces (nSF2) . 
with two closed gapless surfaces (nSFl), and the fully 
gapped superfluid state (SF), respectively. In the case of 
7 = 0, these gapless surfaces collapse into gapless points 
on the /cz-axis, which is consistent with the results of a 
3D Fermi gas with Rashba SOC |16j . 

With the boundaries between different superfluid 
states fixed, we may then map out the typical phase di- 
agrams on the ^ — a plane. As discussed in the previous 
section, the phase diagram reflects the phase structure 
of a trapped gas under LDA. In Fig. |2] we show the 
phase diagrams for pairing states with Q = and under 
an isotropic SOC (7 = 1). As we have anticipated, there 
are two different nodal superfluid states on the phase di- 
agram, with either four (nSF2) or two (nSFl) gapless 
surfaces in momentum space. The stability region of 
the nodal superfluid states increase slightly toward the 
BCS side. Another prominent feature is the existence of 
first-order phase boundaries on the phase diagram. Sim- 
ilar first-order boundaries have been reported for Fermi 
gases in various dimensions and with different forms of 
SOC and effective Zeeman fields [111 US HE] . Apparently, 
these first-order boundaries are due to the competition 
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Figure 3: Phase diagram in the fi — a plane with 7 = 0.5, 
(khtts)'^ ~ —1. The bold solid curves represent the first-order 
boundaries, and the thin dashed curves represent continuous 
phase boundaries. 



between SOC, Zeeman field and pairing. The existence 
of these first-order boundaries implies double-well struc- 
ture in the thermodynamic potential, and requires ex- 
plicit consideration of the phase-separated states for a 
uniform gas. We note that due to numerical difficulties, 
the continuous phase boundary of the normal state is 
determined by the threshold A/h = 10^^. 

To demonstrate the infiuence of anisotropic SOC on 
the pairing states, we also show in Fig. [3] the phase di- 
agram for 7 = 0.5. It appears that the SOC anisotropy 
enhances the stability of the nodal superfiuid states, and 
destabilizes the phase-separated state as it makes the 
first-order boundary shorter. In the case of 7 = 0, we 
would recover the corresponding phase diagram of a 3D 
Fermi gas with Rashba SOC and an axial Zeeman field, 
where there should be considerable stability region for the 
nodal superfluid states. The first-order boundary though 
would persist in the 7 = limit [H] . 



IV. FFLO PAIRING STATES UNDER 3D SOC 
AND ZEEMAN FIELD 

With the understanding of zero center-of-mass momen- 
tum pairing states, we now investigate possible FFLO 
pairing states under 3D SOC and an axial Zeeman field. 
In the weak coupling limit, it is easy to see that the com- 
bination of 3D SOC and the axial Zeeman field makes 
the Fermi surfaces asymmetric in momentum space along 
the z direction. Furthermore, it has been pointed out be- 
fore in different contexts that for systems with SOC and 
asymmetric Fermi surfaces, the conventional zero center- 
of-mass momentum pairing state would become unstable 
against finite center-of-mass momentum FFLO pairing 
states dSl nil [21]. We will show that this is also the 
case here. To see this point more clearly, we first per- 
form a small center-of-mass momentum expansion of 
the thermodynamic potential Eq. (pi around the local 
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Figure 4: Phase diagram for FFLO states in the /i — a plane, 
with (a) (fchtt,)-^ = 0, 7 = 1 and (b) {khas)~^ = -1, 7 = 0.5. 
The bold soHd curves represent the first-order boundaries, and 
the thin dashed curves represent continuous phase boundaries. 



minimum with Q ~ (see Appendix [A] for details) 

n{A,Q,) = no{A) + n,{A)Q, + n2{A)Ql + 0{Ql). (7) 

We find that for finite pairing order parameter A and 
finite axial Zeeman field h, the first-order expansion co- 
efficient Qi is typically finite and is opposite in sign to h. 
Therefore the local minimum in the thermodynamic po- 
tential with Q = is shifted onto the finite plane. For 
h > 0, this suggests that pairing states with zero center- 
of-mass momentum in the presence of 3D SOC are always 
unstable against an FFLO pairing state with < 0. 

We then map out typical phase diagrams on the fi — a 
plane for pairing states with finite center-of-mass momen- 
tum (see Fig. |4]). On the phase diagrams, all the zero 
center-of-mass momentum pairing states are replaced by 
FFLO pairing states with center-of-mass momentum op- 
posite to the direction of the axial Zeeman field {Q^ < 0) , 
which is consistent with our previous analysis. Note that 
for appropriate parameters, pairing states with center-of- 
mass momentum perpendicular to the Zeeman field are 
also stable against pairing states with Q — 0, but are 
metastable against the ground state on the phase dia- 
gram. As we have discussed before, a qualitative under- 
standing here is that the local minima of the thermody- 
namic potential in the Q = case are shifted onto the 
finite Qz plane. When this shift in phase space is not too 
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Figure 5: Typical gapless contours for the FFLO states in 
the kx — kz plane. The parameters are: (a) akh/h = 0.38, 
li/h = 0.46; (b) aku/h = 0.7, ^l/h = 0. For both figures, 
7 = 1, (fchfl,)-! =0. 



large, many properties of the resulting FFLO states are 
similar to those of the original zero center-of-mass mo- 
mentum pairing states. For instance, we still have first- 
order phase boundaries between different FFLO states or 
between FFLO state and the normal state; we still have 
nodal FFLO states with gapless surfaces in momentum 
space (see Fig. [s]); and we have a continuous boundary 
between a gapless FFLO state (FFLO) and a fully gapped 
FFLO state (g-FFLO). These FFLO states, the g-FFLO 
state in particular, are characteristic of pairing states in 
the presence of SOC and Fermi surface deforamtion. For 
example, in the weak coupling limit, one can clearly see 
that the g-FFLO state features an SOC-induced single- 
band pairing, which necessarily has finite center-of-mass 
momentum due to the Fermi surface asymmetry. As 
the Fermi surface asymmetry decreases with increasing 
chemical potential /i and/or SOC strength a, we find 
that the center-of-mass momentum of the g-FFLO state 
decreases in magnitude with increasing /i or a. We expect 
that similar FFLO pairing states should exist in systems 
with general forms of SOC and Fermi surface asymmetry. 

Finally, we emphasize that due to the existence of first- 
order phase boundary on the phase diagram, one must 
explicitly take phase-separated states into account to get 
the correct ground state of a uniform gas. As we have dis- 
cussed previously, in a uniform gas, the total particle den- 
sity is typically fixed by the number equation. However, 
across the first-order phase boundary, the total particle 
density of the ground state does not change continuously 
(see Fig. [6](a) ) . This suggests that for certain total num- 
ber densities, the number equation and the gap equation 
could not be solved simultaneously, unless a phase sep- 
arated state is considered. More explicitly, we show in 
Fig. |6jb) the phase diagram on the n ~ a plane, where 
n — —dQ/d^ is the total number density. It is clear 
that there exists a phase-separated region (PS), in which 
the number equation and the extrema condition of the 
thermodynamic potential, i.e. the gap equation, cannot 
be satisfied simultaneously. This is direct evidence that 
one must explicitly consider the possibility of phase sep- 
aration for a uniform gas. Similar argument should also 
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Figure 6: (a) Total particle number density across the first- 
order piiase boundary for fixed SOC strength akh/h = 0.4. 
(b) Phase boundary on the n — a plane. For both figures, 
7 = 1, (fchOs)"^ = 0. The unit for the particle density no = 



apply for calculations in the canonical ensemble, where 
the total particle number is fixed. 



FFLO states in a polarized Fermi gas. With the recent 
proposals for realizing 3D SOC in ultracold atom gases, 
our work is helpful for the future experimental investi- 
gations with 3D SOC in ultracold Fermi gases, and pro- 
vides valuable information for the general understanding 
of pairing physics in spin-orbit coupled fermionic sys- 
tems. 
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Appendix A: Instability of BCS state against FFLO 
pairing 

In this appendix, we present some details regarding the 
expansion of thermodynamic potential Eq. ([s]) around 
the BCS state with respect to small center-of-mass mo- 
mentum Qz. For the BCS state, we denote the eigen- 
values of the matrix M^~^ in the effective Hamiltonian 
Eq. (3) as E^^{Q — 0), where i^, A = ±. By considering 
a small center-of-mass momentum Q^, the four quasi- 
particle (hole) dispersions can be written as 



0(Q^)(A1) 



V. SUMMARY 

We have studied the pairing states in a 3D Fermi gas 
under 3D SOC and an effective axial Zeeman field. We 
show that the interplay of SOC, Zeeman field and pairing 
superfluidity lead to rich phase structure with exotic pair- 
ing states. In particular, over large parameter regions, 
the ground state of the system is an FFLO state with 
finite center-of-mass momentum. These FFLO states ex- 
hibit interesting excitation properties, ranging from hav- 
ing gapless surfaces in momentum space to being fully 
gapped. We argue that these FFLO states are unique 
to the coexistence of SOC and Fermi surface asymme- 
try, and are qualitatively different from the conventional 



To determine the expansion coefficients (S's, we employ 
the identity 

tr(A") = ^Ar (A2) 

i 

for n < 4, where Xi is the ith eigenvalues of matrix A. 
By matching coefficients of terms involving Qz and 
on both sides the equation above, we can write down 
two sets of linear equations for S^^^^ and k jyi respec- 
tively. After some straightforward algebra, the expansion 
of thermodynamic potential Eq. ([5| with respect to small 
Qz can be obtained via substitution of the resulting dis- 
persions. 
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